Single particle in quantum gravity 
and Braunstein-Ghosh-Severini entropy of a spin network 
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Passermi and Severini have recently shown that the Braunstein-Ghosh-Severini (BGS) entropy Sp = 
— Tr [pr log (Or] of a certain density matrix pr naturally associated to a graph Y, is maximized, among 
all graphs with a fixed number of links and nodes, by regular graphs. We ask if this result can play a 
role in quantum gravity, and be related to the apparent regularity of the physical geometry of space. 
We show that in Loop Quantum Gravity the matrix pj- is precisely the Hamiltonian operator (suitably 
normalized) of a non-relativistic quantum particle interacting with the quantum gravitational field, if 
we restrict elementary area and volume eigenvalues to a fixed value. This operator provides a spectral 
characterization of the physical geometry, and can be interpreted as a state describing the spectral 
information about the geometry available when geometry is measured by its physical interaction 
with matter. It is then temptir\g to interpret its BGS entropy Sf as a genuine physical entropy: we 
discuss the appeal and the difficulties of this interpretation. 



I. INTRODUCTION 

The content of this paper is twofold. First, we con- 
struct a theory describing a single quantum particle in- 
teracting with the quantum gravitational field, in the 
approximation where the particle moves slowly. We 
focus on the quantum dynamics of the particle, and 
study how this is affected by the quantum discreteness 
of spacetime. We use Loop Quantum Gravity (LQG) 
canonical methods llHS]. This simple model might be 
of help in view of the general problem of describing the 
behavior of matter on a quantum spacetime. But it can 
also teach us how the quantum properties of geometry 
can be measured by a material apparatus. Interestingly, 
the Hilbert space of the particle turns out to be natu- 
rally dependent on the geometry, and to be restricted to 
states with support on nodes of the spin network that 
describes the state of spacetime. Our main result is the 
construction of the Hamiltonian operator that governs 
the particle dynamics. The operator depends on the 
state of the quantum gravitational field, since the par- 
ticle energy is both a function of the particle degrees of 
freedom and of the geometry. 

The second part of the paper is based on the obser- 
vation that in the approximation where we can take 
the elementary area and volume quantum numbers 
of the spin-network describing the quantum gravita- 
tional field all equal, the Hamiltonian of the particle, 
suitably normalized, is the Braunstein-Ghosh-Severini 
(BGS) density matrix pi of the graph F of the spin net- 
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work describing the quantum spacetime In the con- 
text of graph theory, Braunstein, Ghosh and Severini de- 
fine a "von-Neumann entropy" Sp = —Tr[pY log Py] as- 
sociated to a graph F. In a recent paper [5], Fasserini and 
Severini show that among all graphs with a given num- 
ber of nodes and links, this entropy is maximized by reg- 
ular graphs: graphs where the links are most uniformly 
distributed among the nodes. This implies that -in the 
approximation given-, the BGS entropy of the (normal- 
ized) Hamiltonian is maximized by spin networks with 
a uniform graph. In other words, there exists a notion of 
"entropy" associated to quantum states of space, which 
is maximized by uniform states. It is natural to ask if this 
graph-theory result can be used in the context of quan- 
tum gravity, and if it has a physical interpretation. It is 
even tempting to try to use the Passerini-Severini theo- 
rem, or some extension of the same, to justify the appar- 
ent uniformity of physical space as it is observed via its 
interaction with matter. We discuss this possibility, and 
its difficulties, in the final part of the paper. In particular, 
we argue that the BGS entropy cannot be directly taken 
as a measure of gravitation entropy, but it may indirectly 
indicate a possible definition of such entropy. 



II. PARTICLE ON A QUANTUM GRAVITATIONAL 
FIELD 



A. General Relativity with a particle 

In this section we construct the quantum theory 
of General Relativity interacting with a single par- 
ticle. In the conventional (ADM) framework, the 
canonical phase space of this theory is coordinatized 
by {qabix), Ti'^ix), X«, P«) where qahix) is 3d metric 
{a,b = 1,2,3) with conjugate momentum tz''^{x) and 
(X", Pa ) are the position and the canonical momentum 
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of the particle, x = {x"} denotes the coordinates of a 
3d constant-time surface. The dynamics is given by a 
Hamiltonian constraint C{x) = Hadm{^) + ^^{x,X)Pq 
where tiyjM is the Hamiltonian constraint of pure grav- 
ity and Pq is the energy of the particle. 

We take the approximation in which the velocity of 
the particle is small^ that is we take ^ and Pq = 



m + 2^ . In this approximation. 



C{x) = HM)M{x)+S^{x,X)m- 



S\x,X)^. (1) 



We focus on the last term, the Hamiltonian of the 
particle on the given gravitational field. Here P^ = 
q"^ [x) PaP}, ; notice the presence of the metric tensor, 
which gives the interaction between the particle and the 
gravitational field. Smeared with the Lapse function 
N{x) the Hamiltonian reads 



H 



— JdxN{x)S^{x,X) q'\x)PaPi, 
q (X) Pa Pb . 



(2) 



In what follows we take N{x) = 1 . 

In view of the quantization it is convenient to rewrite 
the theory in terms of Ashtekar variables. We shift to 
the variables ( E'"'(x), A^(x), X", ) where A[{x) is the 
Ashtekar— Barbero connection [3], and E'"(x) is the cor- 
responding conjugate momentum; this is the densitized 
triad, related to the 3d metric by E"E^^ = qq"^ , with 
/ = 1, 2, 3, where q = det (j^^. In terms of these variables, 
the Hamiltonian of the particle reads 



H 



1 E'"(X)E'"(X) 



Im 



q{X) 



PaPb- 



(3) 



This is the Hamiltonian that we want to define in the 
quantum theory. 



B. Quantum theory 

The state space "H = "HiQsf^i'Hp of the quantum the- 
ory will be the tensor product of the gravity state space 
"Hyj; and a space of particle states T-Lp. In the first, we 
chose a spin network basis that diagonalizes areas and 
volumes. That is, we take the LQG state space 'Hiqg, 
spanned by a basis | s ) labeled by spin network s. A 
spin network s = (T, ]£, t/„ ) is a graph F immersed in IR'', 
colored with quantum numbers on each link £ and 



quantum numbers v„ on each node n. Here v„ are vol- 
ume eigenstates. That is, we choose a basis | v„ ) in the 
space of the intertwines at each node, that diagonalizes 
the volume of the node. These states are linear combina- 
tions of states spanned by the conventional intertwiner 
basis. The states | s ) are orthonormal. 

In the particle state space T-Lp we chose a (possibly 
generalized) basis | x ) that diagonalizes the particle's 
position X. We write states forming a basis in = 
'HiQs^'Hp in the form 



s, X , 



s) (g) I x) en 



(4) 



The definition of the appropriate scalar product on Hp 
requires some care. In usual non-relativistic wave me- 
chanics, the scalar product is usually written in the form 



dx !/;(x)<^(x) 



(5) 



but in this expression the coordinates x are the Cartesian 
coordinates of space, that contain metric information. In 
general coordinates, the above expression reads 



dx^/q xp{x)(p{x) . 



(6) 



Taking this as the definition of the states | x ) is equiva- 
lent to write in "Hp the resolution of the identity 



dx I X ) ( X I 



that is, to have states normalized as 



1 



5{x,y). 



(7) 



(8) 



Equation l|7ll is the appropriate resolution of the identity 
in T-Lp. However, in the theory we are interested in, ^ is 
an operator on ?^iq3. Therefore the scalar product in Hp 
depends on "Hiqg-^ On H = Huj^^Hp we can sandwich 
the right hand side of ^ between two states ( s | and 
I s ) in ?^iQ3. In this way, we have that the operator ^/q 
is replaced by its expectation value in | s ) : 

( s I ydx |x)(x||s) = J dx {s \^/q\ s ) | x ) ( x | . (9) 

But the volume operator vanishes everywhere except at 
the nodes: 



^Jqix) I s) = J2 v„S{x,x„)\s) , (10) 

n6N(s) 



^ This does not imply that we take c co in the theory | 6]. It only 
means that we make measurements in frames where the particle 
moves slowly. 



^ That is, by itself Hp is a vector space but not a Hilbert space. It is 
H = 'Huj:;® 'Hp which is a Hilbert space. 
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where we have indicated as N{s) the set of the nodes 
« of s that have nonvanishing volume eigenvalue v„. 
Using the last equation in the previous one gives 

(s|lp|s) = J2 Vn\x„){x„\, (11) 

«eN(s) 

This gives the resolution of the identity in 'HiQi^'Hp 

t = J2 v„\s,x„){s,x„\. (12) 

^ neN{s) 

It follows that all wavefimctions i/'(x) = {x\ip) that 
have the same values xp{xn) = {x„\ip) on the nodes 
are to be identified as elements of the Hilbert space of 
the theory. The relevant Hilbert space is spanned by the 
states \s,x„) , where n E N{s), satisfying 



^S,Xn\s , X„i ) — ^ss' ^nn' 



(13) 



The quantum states of the particle cannot be considered 
independently from those of the geometry. Remarkably, 
the same conclusion has been reached in [7], in the con- 
text of quantum cosmology, from a rather different set 
of arguments. 

For later purposes, it is also convenient to define the 
states 

\x) := i \x) (14) 

that satisfy 

= J dx q\x){x\ = J2 ^n\^){^ \ ' (15) 



^ S, Xn I S , X^f ^ 



neN{s) 



2 ^ss'^nn' 



(16) 



and the orthonormal states \x) := y^q{x) \ x) . The 
definition of these other bases is only for convenience of 
notation and is not strictly needed; in particular, it is not 
meant to conceal 1/v factors, since inverse volume fac- 
tors are not harmful due to our definition of the Hilbert 
space. This completes the definition of the Hilbert space 
of the theory. 

The quantum theory is obtained by promoting the 
phase space variables to operators on this Hilbert 
spaces. In particular the particle momentum operator 
is Pa = —ih Da , being Da the co variant derivative, and 
the densitized triad operator is defined'^ by [1,^ 



^ In the recent literature, this operator is mainly used as a flux opera- 
tor, that is, smeared with a two surface. As we shall see below, here 
the space index a turns out to be contracted with a space derivative 
rather than with the normal of a surface; hence we need the orig- 
inal form of the operator. The expression tlTl is directly obtained 
from the connection representation acting with £"'(x) = kTi f 

on the holonomy hi{A) = V exp J^dti"{t)Ai,{e{t))T', which is the 
kernel of the integral transform relating the loop and the connection 
representations. 



E'"'(x)|s) =KnJ^ dte''{t)S^{xJ{t))\s,T') (17) 



where £''{t) is the tangent to the link in the point t, the 
sum is over the links £ of the spin network s and | s, T/ ) 
indicates the spin network | s ) with a grasp in the point 
£{t). Here k := S/rGNewton- We also recall for later con- 
venience that I s, t't' ) = + l)\s) if the two grasps 
act on the same point, situated on the link £. 



C. The Hamiltonian operator 

We now study the quantum operator in H corre- 
sponding to the phase space function Q. We chose the 
following symmetric ordering of this operator 



H=— —PaE"{X)E''iX)P, — . 



2m ^ 



Vq 



(18) 



Sandwiching this operator between two resolutions of 
the identity gives 



H 



2m 



dxdy\x){x\ PaE'"{X)E'"{X)Ph \y){y 



2m 



dxdy\x) DaE 



The product of field operators E'"{x)3{x,y)E^\y) at the 
same point is ill defined. This is to be expected and we 
can interpret it as a manifestation of the infinities due 
to the the pointlike nature of the particle. To deal with 
these divergences, we point-split the operators by reg- 
ularizing the delta distribution. We replace it with a 
smearing function /r(x, X) where is zero if |x — X| > 
R and /r = = j^j^ if |x — X| < R. Here R is a small 
real number, representing a regularization (R can be in- 
terpreted as the coordinate size of the particle), and the 
distance | x — X | is the coordinate distance (equivalently: 
the distance in a fiducial background metric rjj^j,, which 
does not affect the final result). This gives the regular- 
ized Hamiltonian 



DaE"{x)fT,{x,y)E^\y)D, 



or 



:i|H|y) 



2m 



q^{x)q^{y) 



DaE''\x)fR[x,y)E^\y)D^ 



2m 



q{x)q{y) 



(19) 



(20) 



Let us now fix a spin network s and study the matrix 
elements of this operator between two states \s,ip) and 
\s,(p) , where ip,cp eHp 

{ip\H,\cp) = {s,ip\H\s,^). (21) 

Inserting two resolutions of the identity (O 

{ip\Hs\cp) = {s,ip\tpHlp\s,(p) (22) 
= J dxdy{s,ip\q\x){x\H\y){y\q\s,(p) 
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Inserting the Hamiltonian | |20|| into this expression gives 



dx dy {s,ip\x] 



2m 

xDaE'''{x)fR{x,y)E''\y)Dj,{y\s,cp] 



— I dxdy daip{x) dbcj){y) fR{x,y) 



x{s\E^-{x)E^\y)\s] 



(23) 



where ip{x) = {x\ip) and we have integrated by parts. 
The integration by parts does not give a boundary con- 
tribution as we can assume the wave function of the 
particle to vanish at infinity (or, more simply, we can 
take the space to be compact.) Inserting the explicit ex- 
pression i flTt for the operators E"' (x) into this expression 
gives 



:^\Hs\cp] 



{xh)^ dxdy daip{x)di,(l){x)fR{x,y) 



- !l 

2m ' 

X ^ /'ds£''(s)^(x,£(s)) 

if.' 

X j^dte\t)S{y,i'{t)){s\s,T'4 



2m 



^ ds dtdsip{£{s))dtcpie'{t)) 



X fR{£{s),e{t)) {s\s,t14 



(24) 



where we have used i''{s)daXp{x) = ds'ip{i{s)). In the 
limit where R is small, the only contribution to this 
sum comes when £ = ^' or when £ and £' meet on a 
node. This second case however gives a lower order 
contribution to the sum. Furthermore, in the same limit 
I s, T^T^, ) = ]({]( + 1) I s ) . Bringing all this together we 
have 



2m 



ds dt 



xdsW{s))dt<p{£{t))fR{£{s),£{t)).{25) 



It is time to deal with our regularization. To this aim, 
recall that the Hilbert structure of the state space of the 
particle depends only on the value of the wave func- 
tion on the nodes. The value of the wave function along 
a link that connects two nodes is physically irrelevant. 
Consistently, the expression /^ds dsXp{£{s)) contained in 
the equation above suggests that this equation should 
only depend on the difference 

Ai^p := ^ds d,ip{£J^) = ip{£j) - ip{£j), (26) 

where £j- and £j are the initial and final points of the 
link £. This is the case for large R, but not for small R, 
where the regularization function fR(£{s), £{t)) cuts the 



integrals off. However, we can exploit the freedom in 
choosing the wave function ip to go to a limit where the 
overall integral does not depend on the regulator. This 
can be done, for each R, by choosing ip{x) to be con- 
stant around each node, so that the region along the link 
where it varies is always smaller than R. In the limit, 
dsip{£{s)) Aiip S{s,Sq), where Sg is an arbitrary point 
on the link. Doing so, we have in the limit 



\xp\Hs\cp) = 



2mVv 



Y^jfiie + l)AiipAe<p. (27) 



We absorb the dimensionless divergent factor Vr — > oo 
into a renormalization of the mass^, defining m* = m Vr, 
and, using also £p^ = hGjsiswton/ we have finally 



2/74 
PI 



2m* 



J2 k{ii + '^) ^e^^e<P- (28) 



Equivalently, if we define the particle states 

\l) ■■= \ej)-\£j) (29) 

we can write the Hamiltonian in the total Hilbert space 
T-L in the form 



H 



{8nhf£^ 
2m* 



E + 1) \s,£){s,£ \ ■ (30) 



This remarkably simple operator is the Hamiltonian op- 
erator of the particle on a gravitational field, and the 
main result of this first part of the paper. We under- 
stand that the derivation above is somewhat acrobatic 
and the regularization a bit brutal. But the final form of 
the operator is simple, convincing, and has the required 
properties, including the correct low energy limit. In the 
appendix, we provide a simpler and more heuristic al- 
ternative derivation of the same operator. 

Let us see some properties of this operator. First, 
the dimension of the first factor is [{8nh)^£p^/2m*] = 
[Energy] x [Length]^; the states | x„ ) have dimension of 

[Length] ^'^ from | [T6l l. Thus H is indeed an energy. 

Second, H is positive semi-definite. This can be 
shown by noticing that it is a linear combination of pro- 
jector operators with positive coefficients. Therefore no 
state can have a negative expectation value. 

Next, the ground state of H corresponds to the case 
where the particle is maximally delocalized. In analogy 



* By this we simply mean that we absorb the infinity into (any) one 
of the constants. An infinite factor is to be expected in this con- 
text, because the particle Hamiltonian we started from is ill defined 
on a discrete geometry due to the product of the triad operators at 
a point. Physically, of course, we can say that it is the limit that 
defines the continuous derivative that must be replaced by a finite 
difference, on a discrete geometry. 
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with the vacuum state of free particle, let us define the 
ground state xpQ [x„ ) = ( Xr, | i/^o ) = Cq constant over the 

I Xn ) states. (That is, its value doesn't depend on n). 

Then ( H i/^o ) =0 and therefore H | ipo ) =0. However, 
if the particle is in this ground state, the probability of 
finding it in different nodes is not the same for all nodes, 
as the states | x„ ) are not normalized. This probability 
is rather 



(31) 



That is, in the ground state the particle is more likely to 
be foimd on nodes with a larger volume. 

It is interesting to notice that if we consider a state en- 
tirely localized on a single node n, the expectation value 
of its energy turns out to be 



; S, Xn |H| S, Xn ) 



2m* 



E^O'^ + l) \ {s,Xn\s,X'n] 

e.en 



(32) 



ten 



If we assume that the volume of the node has Planck 
scale, v„ « then ( s,x~„ \H\ s,x~n ) « 2m*JL^,,„ = h 
i.e. the energy is given by Planck constant over the 
Schwarzschild radius of the particle rg . 

Finally, we observe that the form l|30t of the Hamilto- 
nian operator that we have found is not too surprising. 
As we show in the Appendix, there is a simple minded 
way of obtaining it, starting from the discretization of 
the Schrodinger Laplacian of the particle on a lattice, 
and writing the lattice geometry in terms of areas and 
volumes. 



oriented, differently with respect to the general case in 
LQG.The graph is constructed with the following rules: 
two nodes are connected at most by one link and a link 
cannot close on the same node where it starts. Two 
nodes are said to be adjacent if they are connected by 
a link. One can construct an adjacency matrix A(r) 
whose elements are [A(r)]„^,„ = 1 if {n,m\ e L(r) and 
[A(r)]„ y = 0, otherwise. The number of links attached 
to a node n gives the degree of the node n e V(r), de- 
noted by An- The sum of the degree of each nodes is the 
degree-sum of the graph dp = I]„£N(r) ^n- The degree ma- 
trix A(r) is a diagonal matrix, having on the diagonal 
the dee ree dn correspondent to each node. The combi- 
natorial Laplacian matrix L{T) is a positive semidefinite 
n X n matrix, defined as the difference between A(r) 
and A(r). Once one divides Ly by the degree-sum dp 

the resulting matrix pr := = MHT)) 



is Hermi- 



tian, positive semidefinite and unit-trace, and because 
of these properties it is said to be the density matrix of 
the graph. 

The BGS entropy of the graph T is defined^ as the von 
Neumann entropy of its density matrix: 



S(r) = -Tr [prlogpr]- 



(33) 



In j^, Passerini and Severrni have shown that at fixed 
number of nodes and vertices, the BGS entropy is max- 
imized by regular graphs. Regular means here that the 
distribution of the degrees of the nodes is uniform. 

Does this result have physical interpretation, in the 
context of LQG? 



B. The Hamiltonian operator codes 
the information about the geometry 



III. BGS ENTROPY OF GRAVITY 

As anticipated in the introduction, the Energy oper- 
ator (|30l l that we have found is strictly related to the 
BGS density matrix of the graph of the spin network, 
and its BGS entropy is maximized by regular graphs. In 
this section, we begin by recalling some basic facts of 
BGS graph theory, we illustrate the mathematical rela- 
tion between this theory and LQG and we investigate if 
the BGS results can be interpreted physically. 

A. BGS entropy of a graph 

We begin by reviewing the notion of entropy on a 
graph as introduced in [4]. A graph T = {N,L} is de- 
fined just as a set of nodes N(r) = {1,2, ...,n} and a 
set of couples of nodes L{T), namely links connecting 
nodes. The graphs under consideration are said to be 
simple and undirected: links are not colored and not 



The properties of a physical system are revealed by 
the interactions of the system with its surrounding. The 
geometry of spacetime is a physical system; it interacts 
with the rest of the world by affecting the energy mo- 
mentum tensor of matter, or, in the canonical picture, its 
energy. Therefore the energy of matter can be seen as an 
interaction Hamiltonian between matter and gravity. In 
other words, the only tool we have to detect the space- 
time geometry is via its effect on our material measuring 
devices; this effect is coded in the way the metric tensor 
enters the matter Hamiltonian. 

In the gravity+particle system considered above, a 
measurement of the particle's energy measures particle 
as well as metric properties. In non-general-relativistic 
physics, we view the spacetime geometry as given; then 
the energy is just an observable of the particle. But this 



^ Following the conventions of information theory, we set logarithm 
to be base 2 and log = 0. 
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energy is also an observable of the geometry: we get in- 
formation on the geometry by measuring energy levels 
and energy eigenstates of matter on this geometry. 

Alain Connes has long stressed the relevance of this 
spectral point of view on the geometry of physical space 
ISl]. To put it pictorially: the only information we have 
about the universe is via the spectra of the light we re- 
ceive from the sky; or, even more simply: our eyes 
only see spectra! The information about the geometry 
of space which is available to matter is therefore coded 
precisely into the matter's Hamiltonian operator. Ge- 
ometry is nothing else than what determines the Hamil- 
tonian operator of matter. 

From this perspective, the Hamiltonian operator (|30] l 
constructed above codes the information about the grav- 
itational field that is available to the particle. As an op- 
erator on 'Hp, indeed, Hg := ( s |H| s ) is a fimction of the 
state I s ) of gravity. Different states | s ) determine dif- 
ferent operators. This operator can be reconstructed, or 
"measured" by a sequence of energy measurements on 
known particle states.^ From this perspective, the role of 
operator and state are -in a sense- reversed: each parti- 
cle state defines a measurement that allows us to deter- 
mine ("measure") the operator H, and therefore mea- 
sure the geometry. 

Let us normalize the operator by dividing it by its 
trace, that is, define the imit trace operator 

pc:=Hs/TrHs. (34) 

The expectation value of this operator'' on an arbitrary 
statistical state pp = J dxdx' pp{x, x') \x){x' \ on Tip of 
the particle is 

E = Tr[pG|Op]. (37) 

This is an interesting expression. Its signification is 
clear: it gives the (normalized) mean energy for a cer- 
tain state of particle and gravity. For a given gravita- 
tional state, Pg can be seen as an operator measuring 



* In general, the diagonal matrix elements of an operator H in a ba- 
sis \n) can be simply measured as the expectation value of the 
basis states: H„„ — Tr[H\n){n\]; while the nondiagonal ma- 
trix elements are the difference between the expectation values 
of the quantum superposition and the statistical superposition of 
the corresponding basis states: i-f„„, — Tr [H ^ " ^ "' ^ ( " I + ( "' I ] _ 



In fact, this construction is more general. Consider a general 
statistical state of the geometry, defined by a density matrix p = 
Ess'(/'g)ss' |s)(s' I on'HLQG- This as well determines a Hamiltonian 
operator on 

Hp =TrLQG[Hp], (35) 
(Tr LQG is the trace in 'Hlqg) and a corresponding density matrix on 

|0g :=Hp^/TrHp^. (36) 

The map p ^ pc so defined sends a density matrix in 'Hlqg to a 
density matrix in 'Hp. We get back to the previous case if the p = 
|s>(s|. 



a property of the particle state: its energy. The other 
way around, if we assume that we can independently 
measure the state of the particle, then each given particle 
state p-p can be seen as an operator measuring the "state" 
of the gravitational field coded by pQ. The "density ma- 
trix" Pq can therefore be interpreted as a representation 
of the state of the quantum gravitational field; more pre- 
cisely, of its features that are accessible to the particle via 
spectral measurements. 

Let us now take the simplifying assumption of disre- 
garding volume and area eigenvalues in the definition 
of a spin network.^ For instance, let s be a spin network 
with i( = i and v„ = v. Then Pq determined solely by its 
graph F. And it is easy to see that it is precisely the BGS 
density matrix of the graph^ 

Pg = Pi- (38) 



C. The BGS entropy in gravity 

We have seen above that the BGS density matrix of a 
graph is an expression of a (possibly statistical) state of 
the geometry, as seen by matter. What is then the BGS 
entropy of this state? 

First of all, as shown in j^, this entropy is a quanti- 
tative measure of the regularity of the graph. As such, 
it is certainly a useful physical tool. But can we push 
the analogy with statistical mechanics further, and inter- 
pret the BGS entropy as a genuine physical entropy? At- 
ternpts to define the entropy of the gravitational field are 
in iMEl and llllil, but the problem of defining this 
entropy is completely open; for a discussion, see |10]. 
Can the BGS entropy provide a further tool for dealing 
with gravitational entropy? If we could do so, then the 
Passerini-Severini result that relates this entropy to reg- 
ularity could have a wide reach, and we could try to use 
it in order to argue that the observed uniformity of phys- 
ical space has a statistical origin. The possibility is very 
tempting. 

We have found a structure remarkably similar to that 
of standard statistical mechanics. First, pQ has all the 
mathematical properties of a statistical density matrix. 
Second, it has also the interpretation of coding the quan- 
tum state of gravity, including, possibly, the statistical 
quantum state. We may interpret the trace from T-L to l-Lp 
that defines pc in | |35] l in footnote [7] as a way of tracing 
out gravitational degrees of freedom not seen by mat- 
ter, giving a statistical character to pc even when deter- 
mined by a pure p. Finally, | [37t can be interpreted by 
saying that if the particle state is known pp is an operator 



^ Alternatively, the BGS construction can probably be generalized to 
include spins and volumes. 

' Therefore pp can be written in the form <30t : this provide a straight- 
forward proof that it is positive semi-definite. 
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that measures the statistical state pc- Since this structure 
reproduces that of statistical mechanics, one is tempted 
to say that pc is a genuine statistical density matrix and 
therefore its BGS entropy 

Sbgs = -Tr [lOGlogpc] (39) 

is a physical entropy for the gravitational field. If it looks 
like a duck, swims like a duck and quacks like a duck, 
then it probably is a duck. 

But not necessarily. The problem is that pc is not a 
statistical density matrix. Its eigenvalues are not prob- 
abilities of outcomes. They do not measure numbers of 
microstates. They are (normalized) matter energy levels 
in a certain geometry. Accordingly, we do not see how 
the BGS entropy (|39] l could be interpreted as a physical 
entropy. There are no arguments that we can see, such 
as an ergodicity hypothesis, or Bayesian equiprobabil- 
ity prior, motivating the expectation that we should find 
physical systems maximizing it. 

In quantum mechanics, a pure state | ip ) determines 
an observable, \'^){ip\ whose interpretation is to check 
whether or not the system is in the state \ip) ■ A den- 
sity matrix p = X]/ Pi I ) ( I can also be interpreted 
as an observable, describing the operation of choosing 
with probability to check whether the system is in the 
state \'(pi) ■ The other way around, we can associate to 
any operator A with unit trace a statistical state p^ such 
that on any state the above operation gives the same re- 
sult as the mean value of A. But the analogy stops here, 
because the eigenvalues of A are possible measurement 
outcomes; while the eigenvalues of p are probabilities of 
observing one or another state. In the present case, the 
eigenvalues of pc are still the possible individual out- 
comes of the energy measurement, and have no proba- 
bilistic interpretation that we can see. 

In fact, one can define a density matrix and an entropy 
based on the Laplacian operator H constructed here, 
and therefore indirectly on the Laplacian of a graph, 
which capture the idea of gravitational entropy of the 
distribution determined by measuring the gravitational 
field via the energy of a particle. To this purpose, fix 
the state of the particle tohe \ip) E T-Lf, and define the 
operator 

H' = {xp\H\ip). (40) 

on T-Llqc- Suppose we measure a certain average value 
E of the energy for the particle in this state. This de- 
termines a probability distribution p for the state of the 
gravitational field in "Hlqg- Following O, p can be ob- 
tained as the distribution that maximizes Shannon en- 
tropy, with the constraint given by available measure- 
ments. This gives the density matrix 

p = Z-^ (41) 
which maximizes the entropy 

S = -TrLQG[plogp]. (42) 



Here Z = TrLQcle"''"] and }i is determined by E by 
E = TrLQG[H'|D] = —d{\nZ)/dji. Therefore the relation 
between this gravitational entropy and the Hamiltonian 
operator is 

S = log Z + /^TrLQG [H' . (43) 

A part from the prefactors, and the exponent instead as 
the logarithm, the key point here is that the trace and the 
operator H' are in T^lqg and not in "Hp. 

In particular, we can take xp to be the state where the 
particle is concentrated on a single node n. This give the 
entropy of a spin network s 

p{s) = Z-i e-f'^^ (44) 

where ds is the degree of a single node. A distribu- 
tion that is probably more physically motivated, since 
it avoids singling out a node, can be defined by taking 
instead the average of the degrees over all the nodes of 
the spin network. This corresponds to the entropy de- 
fined by the coarse graining obtained by measuring the 
geometry with a particle concentrated on one node, but 
without knowing which one. Then 

pis) = Z-i e-f"^'^, (45) 

where 

The partition function can be easily calculated [18], and 
gives 

z(H) = E^-'^'= = E (^) ^"''^ = (1+^'-''^)' 

(47) 

This is a very simple possible distribution over spinnet- 
works. Its properties will be explored elsewhere. 

IV. CONCLUSIONS 

In the first part of the paper we have constructed the 
Hamiltonian operator of a quantum particle that moves 
slowly on a quantum gravitational field. The operator, 
given in l l30t has a simple and appealing form. Perhaps 
this construction can be used to explore more exten- 
sively the general form of effective gravity-matter cou- 
plings in quantum gravity. A key result of this part of 
the paper is that the quantum kinematics and the dy- 
namics of the particle cannot be constructed in a way 
disjoint form the geometry: the particle wave function 
turns out to have support on the nodes of the spin net- 
work, and its scalar product depend on the volume, 
namely on the spin network itself. This result (see also 
f?]) can perhaps be relevant for understanding the dy- 
namics of matter in quantum gravity. 
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In the second part of the paper we have observed that 
(under a suitable simplifying assumption) this operator 
is the BGS density matrix of a graph related to the state 
of gravity. We have asked if the corresponding BGS en- 
tropy can be interpreted as a genuine physical entropy, 
related to a probability of finding the state. The struc- 
ture of the theory is extremely similar to that of statis- 
tical mechanics: a trace in the full state space of grav- 
ity determines an operator on a smaller space, which 
has all the mathematical properties of a density matrix 
and codes the spectral information about the geometry 
which is available to matter. 



In spite of the close analogy, this mathematical density 
matrix does not seem interpretable as a physical density 
matrix, since it does not have a precise probabilistic in- 
terpretation. Rather, the density matrix and entropy of 
the gravitational field determined by a matter spectral 
measurement are given by (|45] l and ||43] |, respectively. 
Our analysis is incomplete since it disregards the dy- 
namics, namely the restriction on the physical states im- 
posed by the Hamiltonian constraint; but we do not see 
how this could improve the situation. Therefore, as far 
as we can see, the BGS entropy of a spin network state 
does not appear to have a direct probabilistic interpreta- 
tion. But the analogy remains tantalizing. On the other 
hand, the Passerini-Severini result shows that the BGS 
entropy is a useful measure of the regularity of physical 
space. 



Finally, we have considered a tentative definition of 
a probability distribution on the space of the spin net- 
works. This is determined by maximizing the Shannon 
entropy, conditioned to the measurement of the geome- 
try made by measuring the energy of a particle of known 
position. The properties of the resulting distribution 
will be studied elsewhere. 



In closure, we would like to compare our discus- 
sion with other attempts to define gravitational entropy. 
However, we have found remarkably little on the sub- 
ject. The entropy of matter on a given geometry is 
well understood and discussed everywhere in the liter- 
ature; the entropy to be associated with a black hole is 
of course largely discussed as well; but the entropy of 
the gravitational field itself (in the full nonperturbative 
theory) is a notion that, as far we can see, is almost en- 
tirely ignored in the literature. The only discussion we 
are aware of is the one given by one of us in [10, 11]. See 
these papers for a more detailed discussion. (See also 
ill] and [13, 14] on related subjects.) The statistical me- 
chanics of a generally covariant field theory (which is to 
say: of our true world) is a chapter of physics still to be 
built. 
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APPENDIX 

We give here a naive derivation of the form l i5Dl l 
of the Hamiltonian constraint. The energy operator 
of Schrodinger wave mechanics on flat space is H = 

— ^d^dn and its matrix elements are 

{xp\H\cp) = ^Jdxd^da<p. (48) 

Let us approximate this integral by a Riemann sum over 
the points of a cubic lattice of lattice space L, and ap- 
proximate the derivative by a difference. We have then 



n a 

^ H^n+r,)-'PM lp{x„+s)-(p(Xn) ^^g-^ 

The first sum is over the nodes of the lattice, and the 
second over the three directions. Summing over these 
two variables is the same as summing over the links £ of 
the lattice, obtaining 

{tp\H\cl,)=^^Eil{WJ)-W,)) (cp{if)-cpi^,)). 
This can be rewritten in the form 

iVW'P) -2EU ^ei^— —)[— —) 

where = and Vn = are the area of the square 
separating two lattice cells (of the dual lattice), which 
is cut by the link i and the volume of the cell around 
the node n. Now, let us imagine that instead of being 
on flat space the particle is on a curved space. Then 
the area and volumes A( and Vn will become non trivial 
functions of i and n. Finally, in quantum gravity these 
quantities will have discrete eigenvalues fl6l,[T^ . giving 
directly l l30t . 
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